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Abstract. We present here an elementary example, for every 
fixed positive integer k, of a strictly stationary nongaussian sto- 
chastic process in discrete time, all of whose /c-marginals are gauss- 
ian. 
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1. Introduction 

It is well known that, for every positive integer n > 1, there exists 
a probability distribution in R n which is not gaussian but has all its 
(n — l)-dimensional marginal distributions gaussian. (See, for exam- 
ple, Section 10.3 in Stoyanov [2].) Using the finer theory of pathwise 
stochastic integrals and martingale methods, Foller, Wu and Yor [T] 
have shown that, for every positive integer k, there exists a variety 
of nongaussian stochastic processes with continuous trajectories in the 
interval [0, oo) which have the same A;- dimensional marginals as the 
standard brownian motion process. Here we present an elementary ex- 
ample, for every fixed positive integer k > 1, a discrete time stationary 
stochastic process which is not gaussian but has all its (A; — 1) -marginals 
gaussian. However, we do not know how to construct such processes in 
continuous time. 

2. The basic construction 

Let k > 1 be any fixed positive integer and let v be a probability 
distribution in M fe , which is not gaussian but has all its [k— l)-marginals 
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gaussian with mean and covariance matrix identity. For example, we 
may choose v to have the probability density function 



where x = {x\, X2, ■ ■ ■ , x k ). If {X\, X2, . . . , Xk) is an IR fc -valued random 
variable with distribution v, then the sequence X±, X2, . . . , Xi, . . . , X k 
with the i-ih term omitted consists of i.i.d. N(0, 1) random variables, 
for each i. 

Now consider a bilateral sequence {(Xni, X n 2, . . . X nk j, — 00 < n < 
00} of i.i.d Revalued random variables with the common distribution 
v as described in the preceding paragraph. Define 

Y n = X n k + X n+ i + X n+ 2 k-2 + • • • + X n+ k~i 1, —00 < n < 00. 

It is to be noted that the sum of the two suffixes in each summand on 
the right hand side is equal to n + k. 

Theorem 2.1. The sequence {Y n , —00 < n < 00} is strictly stationary, 
(k— l)-step independent with every (k— 1)- dimensional marginal being 
gaussian with mean and covariance matrix kl, I being the identity 
matrix of order k — 1. In particular, {Y n } is ergodic. 

Proof. Fix an integer m and consider the two sets {Y n , n < m} and 

{Y n , n > m + k}. Since Y rn = X rnk + X m+1 fc _H h X m+k ^i 1 and 

Y m+k = Y m+k k + X m+k+ i fc-i H h X m+2k -i 1 and the first suffix in 

the last summand in the definition of Y m is less than the first suffix in 
the first summand in the definition of Y m+k it follows that the two sets 
{Y n , n < m} and {Y n , n > m + k} are independent. In other words 
{Y n } is a (k — l)-step independent process. 
We now look at the column vector- valued random variable 



V<x) = (2tt)-§ { 



1 + Xix 2 ■ ■ ■ x k e 2 
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X n+ i k + X n+ 2 k-1 + • • • + X n+k 1 
X n+ 2 k + X n+ 3 k -l + • • • + X n+k+ i 1 



Y 



n+rn . 



x., 
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and express it as Si + S 2 + S3 where 



Si 



X n+ i k 




+ 



X n +2 k-1 
X n +2 k 





X n +k 1 

X n +k 2 
X n +k k 
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+ ••• + 
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X n +rn+2 1 
X n +m+2 2 



X n +k-l 2 
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+ ••• + 



X 
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+ ••• + 



X 
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In each column on the right hand side of Si or S3 there are at most 
k — l nonzero entries whereas in each column on the right hand side 
of S 2 there are exactly k entries. All the column vectors appearing 
in Si,S 2 ,S 3 together are mutually independent. By the choice of the 
measure u, Si and S3 are gaussian random vectors. Denote by j]) 
the (j — i + l)-dimensional standard normal distribution imbedded 
in M m so that the first i — 1 and the last m — j coordinates are 
when 1 < i < j < m. Similarly, denote by v{\j, k + j — 1]) the k- 
dimensional distribution v imbedded in M m with the first j — 1 and 



4 



K. R. PARTHASARATHY 



the last m — k — j + 1 coordinates for 1 < j < m — k + 1 , assuming 
m > k. Then it follows that the random variables Y n+1 , Y n+2 , . . . , Y n+m 
expressed as a single column vector has the m-dimensional distribution 
v m (in R m ) given by 

v m = M[MD*M[1,2])*...*M[1,A;-1]) 

k]) * u([2, k + 1]) * . . . * v([m - k + 1, m) 
*fi([m — k + 2, m]) * fi([m — k + 3, m]) * ... * fJ>([m, m]), 

for every m > k. Since v m is independent of n it follows that {Y n } is 
a strictly stationary process. Since is nongaussian it is clear 

that v m is not gaussian for every m> k. 

We now observe that Y n , being a sum of k independent N(0, 1) ran- 
dom variables, is an N(0, k) variable with mean and variance k. Now 
consider the pair (Y , Y rn ). If m > k we have already seen that Y and 
Y m are independent. If m < k, we write 



" ^0 " 




Y 

1 m 





Xok + Xi k-l + • • • + X m _i k-m+l 







+ 



y 

-^■m k—m 
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+ 
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Y 

-^m+l k—m—1 
X m+ l k-l 
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+ + 







Xkrn + X. 



k+1 m-2 



Xk-i i 

Xk-1 m+1 
k+m-1 1 



Now the special choice of v implies that Y and Y rn are independent 
^(0, k) random variables. Stationarity of the process {Y n } implies that 
Y ni and Y n2 are independent N(0, k) random variables for any n 1? n2. 
Now consider, for any n\ < n 2 < • • • < rik-i the random vector 



r y i 




Y 




Y 





X ni k + X ni+ i k-l 



X n2 k + x,. 



x„ 



+ x n 



n 2 +l k-l 



-1 k-l 



X ni+ k-i i 

X n2+ k-l 1 



+ ■ ■ ■ + x, 



nfc-i+fc-1 1 



The right hand side can be expressed as a sum of column vectors in 
which the entries in each column are either or an X rs where the first 
suffix r is fixed and the second suffix takes at most k — 1 values from the 
set {1,2, ... ,k}. The different column vectors are independent and by 
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the choice of v each column has a multivariate gaussian distribution. 
Thus Y is gaussian. Since any two and Yj are independent where k > 
2, it follows that Y ni , Y n2 , . . . , Y nh _ x are i.i.d N(0, k) random variables. 
This completes the proof. □ 

Remark 2.2. From the proof of Theorem 12. II it is clear that any k — 1 
of the random variables {Y n } are i.i.d N(0,k). This motivates the 
introduction of the following notion of limited exchangeability. We 
say that a stationary random process {Z n ,—oo < n < oo} is k- 
exchangeable if any Z ni , Z n2 , . . . , Z nk has the same distribution for any 
A;-point set {ni, n 2 , . . . , nk} C Z. The probability measures of all such 
^-exchangeable stationary processes constitute a convex set. One won- 
ders what are the extreme points of this convex set. 
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